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Abstract. We present a dynamical analysis of a classical billiard chain — a channel 
with parallel semi-circular walls, which can serve as a model for a bended optical fiber. 
An interesting feature of this model is the fact that the phase space separates into 
two disjoint invariant components corresponding to the left and right uni-directional 
motions. Dynamics is decomposed into the jump map — a Poincare map between the 
two ends of a basic cell, and the time function — traveling time across a basic cell of a 
point on a surface of section. The jump map has a mixed phase space where the relative 
sizes of the regular and chaotic components depend on the width of the channel. For 
a suitable value of this parameter we can have almost fully chaotic phase space. We 
have studied numerically the Lyapunov exponents, time auto-correlation functions and 
diffusion of particles along the chain. As a result of a singularity of the time function 
we obtain marginally-normal diffusion after we subtract the average drift. The last 
result is also supported by some analytical arguments. 
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1. Introduction: uni-directional billiard channels 

The discussion of classical and quantum dynamics of spatially extended billiard chains, 
either with periodicity or disorder, is a promising field of research with a variety of direct 
applications, e.g. in nanophysics, fiber optics, electromagnetic cavities, etc. It is fair to 
say that studies in spatially extended billiard systems have been under represented as 
compared to a vast amount of work which has been dedicated to billiards on bounded 
domains. Nevertheless, one has to mention several basic results in this type of systems. 
First, one can study the escape rates from finite portions of an infinite billiard chain, 
like the Lorentz channel pQ. Second, one can study classical transport properties, such 
as diffusion and transport of heat along the billiard chains in order to understand the 
dynamical (microscopic) origin of the macroscopic transport laws [21 HUE]. Third, one can 
study the relation between deterministic diffusion in a classical billiard chain, Anderson- 
like dynamical localization in the corresponding quantum chain, and the nature of its 
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spectral fluctuations EH E| • Ariel fourth, there is an interesting effect of localization 
transition in the presence of correlated disorder, which has been studied in the case of 
billiard chains both theoretically [B] and experimentally 
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Figure 1. A schematic picture of the serpent billiard model (a) with one basic cell put 
in a rectangular frame, and index n labeling the consecutive cells. A typical trajectory 
is indicated in (b). 

In this paper we discuss a class of classical billiard channels with an unusual and 
distinct dynamical property, namely uni-directionality of the ray motion along the chain. 
Specifically we focus our study on a particular billiard chain — a channel with parallel 
semi-circular walls which we name as serpent billiard. The billiard under discussion is 
built as a periodic composition of semicircular rings of radii q G [0, 1) and 1, for the inner 
and outer circular arc, respectively. The geometry of this billiard chain and illustration 
of the ray dynamics are shown in fig. ^ By construction the phase space separates into 
two disjoint invariant components corresponding to the left and right uni-directional 
motions, corresponding to two different signs of the angular momentum as defined with 
respect to the origin (center) of the ring of the current billiard cell. Within each cell the 
angular momentum is conserved. Further, it is obvious that upon the transition between 
one cell to another the sign of angular momentum, as calculated with respect two the 
centers of adjacent cells, remains unchanged. Therefore, particles traveling from left to 
right initially will do so forever and will thus never be able to change the direction of 
travel, and similarly for the motion in the opposite direction, so that these two motions 
constitute two disjoint invariant halves of the phase space. Nonetheless, as we shall show 
below, the dynamics inside each invariant half of the phase space may be (practically) 
totally chaotic and ergodic. 

We note that this property of uni-directionality can be proven for a more general 
billiard channel which is bounded by an arbitrary pair of parallel smooth curves. Namely, 
it is easy to prove the following observation. 

Let the billiard motion in R 2 be bounded by two smooth C l curves Cj, j = 1,2, with 
natural parameterizations s —>■ rj(s). The curves C\ and C2 should never intersect and 
they should be parallel in the following sense: for any s G M, a line C(s) intersecting 
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C\ perpendicularly at r\{s) should also intersect C 2 perpendicularly, say at point r 2 (r) 
defining a map r = o~(s). The function a : K — > R should be monotonously increasing 
invertible function, i.e. cr'(s) > for all s, or in other words, the lines C(s) should not 
intersect each other inside the billiard region. 

Then the billiard motion is uni- directional, i.e. the sign of the tangential velocity 
component v-(d/ds)fj(s) stays constant for all collision points of an arbitrary trajectory. 

To prove this observation, it is sufficient to consider two subsequent collisions of a 
segment of trajectory with a velocity of unit length \v\ = 1. We may assume the first 
collision to take place at r\(s) G C\ and write sina := v • (d/ds)r*i(s) > 0. Then we 
consider two possible cases: 

(a) The next collision happens on the curve C 2 , say at the point f^r). Then the angle 
of incidence again writes as sin (3 := v ■ (d/dr)r*2(r). a and (3 are the angles between the 
segment of the trajectory and lines C(s) and £(er _1 (r)), respectively. Since the latter 
two do not cross inside the billiard region, it follows that the sign of a and (3 should be 
the same (positive). We may have (3 = only if a = 0, i.e. when the motion takes place 
along C(s) which is a periodic orbit. See fig. Et- 

(b) Another possibility is that the next collision happens with the same curve, i.e. at 
r*i(s'). Writing sin/3 = v ■ (d/ds)r 1 (s / ) we again observe that the sign of the angles a 
and (3 should be the same (positive), considering that the other ends of the lines C(s) 
and jC(s') at points r 2 (o"(s)) and r 2 (a(s')) should lie on the same side of the trajectory 
segment since the latter should not cross C 2 . See fig. Eb- Thus we have proven that 
C(s) is a family of marginally stable periodic orbits, of vanishing overall measure, which 
separates the phase space of the billiard into two halves of unidirectional motions. 




(a) (b) 

Figure 2. The schematic diagram shows two possible cases (a,b) of subsequent 
collisions needed for the proof of unidirectionality. 

It is perhaps worth stressing that the conditions of parallelism as expressed in 
the statement imply also that the width of the channel should be constant, \ri(s) — 
^2(o~(s)) \ = const. 

We should note that detailed understanding of the dynamics of such a class of billiards 
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may have useful applications, in particular in fiber optics, electromagnetic waveguide 
propagation, etc. In the following sections we shall concentrate on the dynamics of the 
specific serpent billiard model, which we shall analyze in terms of a special version of the 
Poincare map, the so-called jump map. Then we shall describe analytically and numer- 
ically the average transport velocity, deterministic diffusion and correlation functions of 
the model. 

2. Dynamics of the serpent billiards 

Let us consider the Hamiltonian dynamics of a particle in a serpent billiard channel. 
We are considering a classical point particle with a fixed velocity of unit size. Due 
to uni-directionality of the motion, as shown above, we may freely choose to consider 
only forward propagation — in the positive direction of it is shown in the 

example of fig. ^p. The forward dynamics of the billiard can be written in terms of 
dynamics within a given basic cell and a transition to an adjacent basic cell. In order to 
fully describe the dynamics we only need to know a Poincare-like transformation which 
maps coordinates of an entry into a cell to the coordinates of an exit and the time 
spent between entry and exit (the entry into an adjacent cell). Thus we formulate the 
dynamics of our billiard chain in terms of a jump map model ^H]- We shall define the 
jump map more precisely below. 

Let the particle enter the basic cell on the left end at the initial position x and 
with the horizontal velocity component v x and travel to the right end where it exits at 
the position x' and velocity v' x . To clarify the notation we introduce a phase space entry 
section <Sl and an exit section 5r 



S L = {(x,v x ) : x e [-1, -q],v x G [-1, 1]}, <Sr = {{x,v x ) : x 6 [q, l],v x G [-1, 1]}. 



The dynamical mapping of an entry point to an exit point shall be denoted by 



The map G can be expressed analytically since the billiard in the circular ring is 
integrable. Since the lengthy expression is not very illuminating we give its explicit 
form in the Appendix. In order to be able to apply G again, we have to transform the 
current exit position x' to the entry position of the next basic cell by a map S : 5r — > <Sl 



With this transformation the conservation of angular momentum around the center 
of the current cell is broken which implies non-integrability of the model. We should 
mention that our serpent billiard falls into the category of semi- separable systems [16]. 
The propagation of a particle from one basic cell to another can then be stated in terms 
of a single map F : S L — > 



G : <S>l —> <Sr; 




(1) 




F = SoG. 



(2) 
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We will refer to F as a jump or Poincare map and to the phase space Sl as a surface of 
section (SOS). In terms of the map F the dynamics over the whole channel is decomposed 
into spatially equidistant snapshots. In order to maintain the whole physical information 
about the dynamics we have to introduce the time function T : <Sl — > M + , i.e. T(x) 
measures the time needed by a particle to travel from the entry point x = (x, v x ) G <Sl 
to the other end of the basic cell. The pair (F, T) now represents the jump model 
corresponding to our billiard channel. Again, the time function T(x) could be written 
out explicitly, though with a cumbersome expression, so we put it in the Appendix. 
However, we should note that numerical routines for computing the map F and the 
function T are very elementary and efficient. 

As a useful illustration of the gross dynamical features of the model we plot in fig. |3] 
the phase portraits of the Poincare-jump map F for different values of the parameter 
q. We observe that the jump map has in general a mixed phase space with chaotic 
and regular regions coexisting on SOS. We also observe that the phase portraits are 
symmetric in x around the mean radius (l+g)/2. We note that the chaotic component is 
always dominant in size over the regular components and for certain regions of parameter 
q values the regular components are practically negligible so the jump map appears to 
become (almost) fully chaotic and ergodic. One nice example where we were unable to 
locate a single regular island has the parameter value q = 0.6. In fig. 01 we plotted the 
relative area of regular SOS components as a function of parameter q. 

In order to quantify the exponential instability of trajectories inside the chaotic 
component of SOS we have measured the average Lyapunov exponent A (as described 
e.g. in JU]). The result for A as a function of q is shown in fig. We observe that 
the chaoticity as measured by A, being equal to the dynamical Kolmogorov-Sinai 
entropy, is always positive and is increasing monotonically with q. This trend was 
somehow intuitively expected, as the number of collisions within the jump increases 
with q. Namely, the bigger number of collisions between entry and exit sections implies 
less correlation between the angular momenta of the adjacent cells meaning stronger 
integrability breaking. 

Important fingerprints of dynamics, in many ways complementary to Lyapunov 
exponents, are the time-correlation functions. These reflect the mixing property of the 
system which implies decay to equilibrium of an arbitrary initial phase space measure. 
Time correlation functions are also directly related to transport, which is studied in the 
next section, through a linear response formalism. We discuss discrete time correlation 
function of the jump map between two observables 0(x) and which is defined as 



The correlation functions are normalized such that always, C^(0, x) = 1, even for 
observables which are not in L 2 {Si). In the sequel we consider autocorrelation functions 
of very regular observables such as the phase space coordinates, namely C X)X and C VxtVx , 
and the autocorrelation function of the time function Ct,t which is even more interesting 




(3) 
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Figure 3. Gallery of phase space portraits of the jump map F for different values of the 
parameter q. Horizontal axis: x, vertical axis: v x . Each diagram shows 10 4 successive 
iterations of 400 different initial points distributed randomly over SOS. 
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Figure 4. Numerical estimation of the total relative area of regular SOS components 
1 — p, where p is the relative area of the chaotic component of SOS, for the jump map 
is shown as a function of the parameter q. The ratio is obtained by sampling of 1000 
random trajectories starting inside the chaotic component over a phase space grid of 
a size 1000 x 1000. The length of trajectories in number of jumps was 10 6 . 
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Figure 5. Lyapunov exponent A of the jump map F averaged over the dominating 
chaotic component of SOS as a function of the parameter q. The average was computed 
from 1000 trajectories of length 10 5 . 



for two reasons: (i) Ct,t is directly related to particle transport as described in the 
next section, and (ii) T is not in L 2 (5l) as discussed below. Numerical data presented 
in fig. |U1 strongly suggest that correlation function Ct,t typically exhibits exponential 
decay ~ exp(— constr) for most of values of q, except that for small parameter values 
q < 0.3 the initial exponential-like decay seems to turn into an asymptotic power law 
decay ~ ^-const Q n other hand the correlation decay of non-singular observables, 
like C x x and C VxjVx seems to behave as a power law for all values of q. It is interesting to 
observe that the qualitative nature of correlation decay seems to be quite different for 
different classes of observables such as x compared to T . However, in all the cases time 
correlation functions strongly decay which is a firm indication of the mixing property 
of the serpent billiard on the chaotic component. 

The time function T(x) is expected to have a singularity for v x — ► ±1 as it may take 
arbitrary long to traverse the semicircular ring with sufficiently small value of angular 
momentum. It is straightforward to show that this is a square-root singularity 



An example of the structure of the time function for q = 0.6 is shown as a density plot 
in fig. [7^. Another quantity which can illustrate the dynamical behaviour of observable 
T is the probability distribution P(T) of times T(x) for a very long chaotic trajectory. 
Assuming that the system is ergodic on the full SOS, the probability distribution P(T) 
can be written in terms of the constant invariant measure on SOS 




v 



x 



(4) 




(5) 



Singularity (@J implies the asymptotic form of the distribution 

p(t) ~ r 3 



(6) 
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Figure 6. The average auto-correlation functions: (Ct,t(x, t))q (full curve), 
(C xx (x., r))n (dashed), and (C VxVx (x, r))n (dots) as functions of the number of jumps 
t. The calculation is performed using trajectories of length 5 • 10 5 and averaged over 10 4 
initial conditions x over the chaotic SOS component. In the plot for q = 0.6 and q = 0.9 
we insert the exponentials exp(— 0.0508193 r) and exp(— 0.00234154 r), respectively, in 
order ti guide the eye. The noisy plateaus indicate the level of statistical fluctuation. 

for large t. This asymptotic property does not essentially depend on the full ergodicity of 
the map, as we find the same asymptotic behaviour by numerical simulation of P(t) for 
different values of q as shown in fig. [TJd. It is obvious that the only important condition 
for the universal decay of P(t) is that the chaotic component should extend to the lines 
of singularity v x = ±1. The distribution P(t) is very important, because it directly 
connects to the particle transport properties of the channel that are discussed in the 
next section. 

3. Transport properties 

Here we would like to examine the transport properties along our billiard chain in the 
context of the jump model. The basic cells are labeled with nonnegative integer n G Z + 
starting with n = and counting forward to the right. Transported length is measured 
as the number of traversed basic cells. This means that our length of travel is an integer 
but we should keep in mind that the physical length of one basic cell in the direction 
along the chain is 1 + q. 

Let us prepare an ensemble of particles on the SOS Sl of the zeroth basic cell, 
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Figure 7. (a) Density plot of the time function in logarithmic scale at q = 0.6 and (b) 
the distribution P(t) of the values of the time function for different q indicated in the 
figure. A thicker full straight line indicates t~ 3 decay. 



n = 0. The phase space is generally mixed and let the dominating chaotic component 
be denoted by Q C Sl- We observe that each invariant phase space component may have 
its own transport properties and to obtain a clear picture of transport we have to test 
each component separately The transport on regular components (islands) of SOS is 
obviously ballistic since the islands are transporting with a constant and sharply defined 
average velocity. Thus we concentrate on the more nontrivial case of transport on the 
chaotic component, namely in the following we choose a uniform initial distribution of 
the particles over the leading chaotic component Q: /i(d 2 x) = d 2 x/vol(fi). 

The transport of an ensemble of particles is described using the probability 
distribution P(n,t) of particles over the basic cells n G Z + as a function of time t. 
Let us express this distribution in terms of a jump map F(x) and the time function 
T(x). The time spent by a particle to traverse m basic cells starting from the initial 
position x G Sl is calculated as 

TO— 1 

T m (x) = 5>(F«(x)). (7) 

fc=0 

The probability distribution of a single particle with initial coordinate x over the basic 
cells (labeled by n G Z + ) at time t can be written straightforwardly as 

P(n,t,x)= f {5(T-T n (x))-5(T-T n+1 (x))}dT. (8) 
Jo 

By averaging P(n,t,x) over an initial ensemble of particles, defined by 

(/(x)>„ := f /i(d 2 x)/(x), (9) 
Jn 

we obtain the distribution of particles over the cells 

P(n,t) = (P(n,t,x)) n = [ {Pn(r)-p n+ i(r)}dT, (10) 

Jo 
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where we are using the ensemble average distribution of times needed to traverse n basic 
cells, denoted asp n (t), 

p n (t) = (S{r-T n {x))) n . (11) 

Let («)p denote an average over the spatial distribution P(n,t) 

oo 

(/(n))p = £/(n)P(M). (12) 

n=0 

Our aim is to obtain the time-asymptotic (t — > oo) form of P(n, £) where only 
contributions of far lying cells n ^> 1 are relevant. Since the auto-correlation functions 
in our model are decaying very fast, in particular the relevant (Ct,t(x, £))n (see the 
end of previous section), we may employ the central limit theorem and approximate the 
distribution of T n , p n (t), in terms of distribution of T, pi(t) = (5(t — T(x)))q: 

p n {t) = ( pi*.-.*pi )(0, n > 1. (13) 

This approximation is very useful, because we can obtain the whole distribution P(n,t) 
using a single function p\(t) = P(t) that can be easily measured and is already plotted 
in fig. [7b. In the limit n->oowe can treat the variable T n as infinitely divisible [TTj and 
the parameter n as a continuous variable. Then we can approximate the finite difference 
in n in eq. |T0|) in terms of a derivative d/d n and calculate P(n,t) as 

P(n, [ Pn (r)dr. (14) 



dn J 

Basic properties of the transport shall be described by the mean traversed length (n) p 
and the spatial spread of the initial ensemble o\ = (n 2 )p — (n)p. The mean (n)p can 
be asymptotically, for t — > oo, exactly expressed by the formula 

(«>p= / $>„,(r)dT. (15) 
^0 n=0 

From the central limit theorem we immediately obtain the mean velocity as the inverse 
mean traverse time 

(n)p If 00 ,, 
hm t mean = - = / t Pl (t)dt (16) 



t^oo t V jq 

The average t mean , and the minimal time t m i n to traverse a basic cell, as a function of q 
are plotted in fig- El The mean time t mea n is almost linearly increasing with increasing the 
parameter q. This reflects the obvious fact that the travel becomes slower by narrowing 
the channel. The fact that the linear growth of (n) p is indeed given by velocity v is also 
demonstrated numerically for q = 0.6 in fig. EK- It is a little more tedious to obtain an 
analytical approximation for the average spreading width cr^it) = (n 2 ) P — (n)p. 

Essentially we need to control the growth of second moment of distribution P(n, t). 
This distribution is given by eq. ()14)) in terms of Pn (t) which may be asymptotically 
expressed as a convolution of independent distributions P i(t) (fT3*|) assuming suflicienly 
fast decay of temporal correlations as established numerically. As a consequence, p n {t) 
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Figure 8. Minimal time i m i n and average time t mean for traversing a basic cell over 
the chaotic SOS component are shown as functions of parameter g.The curve that 
fits imean is the linear function 2.45863 + 2.48696g. The averaging and minimizing is 
performed on a trajectory of length 10 7 . 



inherits cubic singularity of pi(t), namely p n {t) oc t~ 3 . This heuristic argument would 
suggest marginally normal diffusion er^(t) oc tlogt, not essentially connected to the 
strength of correlation decay but simply as a consequence of singularity of the time 
function. 

This observation can be formalized with a brief calculation. We stress that for large 
times t —>■ oo only cells with labels n ~ v t ^> 1 contribute appreciably to the probability 
distribution P(n, t). In this regime we approximate p n {t) using hence we can easily 
write its Fourier transform as 

/oo 
exp(ite)/(t)dt. (17) 
-oo 

The asymptotics of pi(k) for k — > will determine the asymptotics for p n (t) for long 
times t — ► oo. Thus we write the local expansion of pi(k) around k = explicitly taking 
into account the known asymptotics in time domain, pi(t —* oo) ~ t~ 3 , namely 

logpi(fc) = it mean k + alk 2 (a + \ogk) + 0(k 3 ), k> 0, (18) 

where <Tq and a are positive constants depending on the details of pi(t). Using eqs. 
fll8pi7j) and applying the inverse Fourier transform we find that the limiting distribution, 
namely p n for large n, is given by the formula 

Pn(t) = —g(- — 7 ^ mean ) , orf = alnlogn, for n > 1. (19) 

where g(x) = ^= exp(— x 2 /2) is the standard Gaussian function. By inserting the above 
expression (fTTJ|) into approximation of P(n, t) (fT^j) and expanding in terms of parameter 
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n — vt we get the result 

P(n,t) = —g , al = alv 3 t\og(vt), for n > 1. (20) 

From this analysis we predict that the diffusion of initial ensemble of particles will 
distribute over the basic cells with dispersion growing as tlogt. Thus we have shown 
analytically that our serpent billiard exhibits marginally normal diffusion with a drift. 

We have tested our analytical results by performing extensive numerical 
simulations. An example of of cr^(i) for q = 0.6 is shown in fig. ^jp. We stress that 
our numerical data are indeed consistent with the marginally normal diffusion. 




10 2 10 3 10 4 10 5 10 6 10 7 10 8 10 9 



t 




t 



Figure 9. The average position (n) p (a) and the spread of particles ofj = (n 2 )p — {n)p 
(b) as a function of time t at q = 0.6. The average transport velocity v = l/i m ean, 
where t mcan = 3.94538 is obtained from the data of fig. |H1 The error bars represent 
the standard deviation of the measurement. The statistics is performed over an 
ensemble of 25000 initial points. The fitted line in the inset is given by cr^/i = 
-0.044374+ 5.2636 log 10 t 
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4. Summary and discussion 

In this paper we have analyzed a simple billiard chain, the so-called serpent billiard, 
with a special dynamical property of strictly uni-directional classical motion. We have 
also proven unidirectionality of motion for a more general class of billiard channels with 
parallel walls. The dynamics along the serpent billiard channel is described in terms of 
a variant of a jump model namely the jump-Poincare map between the surfaces 
of section of two adjacent basic cells of the billiard, and the time function i.e. the time 
needed to traverse the basic cell as a function of the position on the surface of section. 
We have shown that the jump map is chaotic with generally mixed phase space, where 
the relative size of the largest chaotic component is generally increasing with decreasing 
of the channel's width. The latter dependence is not strictly monotonic, because of 
bifurcations of regular components, but for narrow channels the chaotic component 
is typically largely dominant. For a considerable range of the parameter (denoting the 
channel's width) the jump map is even practically fully chaotic, ergodic, as no detectable 
islands of regular motion have been found. This does not mean that the islands of 
stability can not exist for typical values of the parameter. We only wish to stress that 
it is easy to find parameter values (such as the case q = 0.6 studied in the paper) for 
which all possible islands of stability are undetectably small for numerical (experimental) 
purposes. Numerically measured maximal Lyapunov exponent shows that the chaoticity 
on a chaotic component is monotonically increasing with narrowing the channel. 

The transport of particles along the channel measured in the number of traversed 
basic cells exhibits a marginally-normal diffusion cr^ ~ tlogt (when the drift term 
is subtracted) due to square-root singularity of the time function. This singularity is 
a consequence of parallel walls, or saying in dynamical terms, it is due to a family of 
marginally stable bouncing ball trajectories bouncing perpendicularly between the walls. 

Besides its interesting and rather exotic dynamical properties, the model and its 
generalizations may also be relevant for real world problems of transport, such as in 
optical fibers or wave-guides. These results may open even more interesting questions 
on the properties of quantum or wave transport of classically unidirectional billiard 
channels. This is a subject of a subsequent publication [TTj . 
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Appendix: Explicit jump map and time function 

Let the particle enter the basic cell at point x = (x,cos0) G <Sl end exit at point 
x' = (x',cos0') G Sr. Here we explicitly write the map G : x — > x', eq. (0), and the 
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time function T(x). Due to conservation of the angulur momentum within a fixed cell 
we have the relation 

T = —x sin = x sin 0' 

so we just have to give an explicit formula for the map g : — > 0' and the time 
function i(0). Then the remaining relation to specify the full map G simply reads 
x' = — xsin0/ sing(0). Let us write some auxilary variables 

(3 = arcsinr, 7 = arcsin(r/g), 

/ r n — a 

a = a + no, n = — - — , 

5 

where |_^J is the largest integer not larger than x, and a and 5 are determined for each 
case separately below. 



If T > q the particle is only hitting the outside wall. Then, writing a = ir — — (3, 

5 = it — 2(3, we have 

0' = (n + 1)5 - 0, 

t = VI + x 2 + 2x cos a + \/ 1 + x' 2 + 2x' cos a' + 2n cos (3 

In the opposite case where T < q, writing 5 = 7 — (3 and At = a/1 + q 2 — 2qcos5, we 
have to discuss two cases: (i) when the particle first hits the inner wall, < 7r/2: 

(n + 1)5 — n odd 
27 — 7r + n5 — n even 



a/1 + x' 2 + 2rr' cos a' n odd 



* = V Q + ^ + 2xgcoso; + nAt + < r - ^ , 

I v 9 + x + 2^ 9 cos a' n even 

and (ii) when the particle first hits the outer wall, > n/2: 

, (n + 1)5 — n odd 
1 7r — 2/3 + n5 — n even 



n o — ^ a, ! y q 2 + x' + 2x' q cos a' n odd 

t = VI + x 2 + 2x cos a + nAt + < V 9 

I v 1 + x' + 2x' cos a' n even 

where a is equal to 7 — in case (i) and n — — (3 in case (ii). 
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